arXiv: 1509.018lOvl [math.AP] 6 Sep 2015 


Stable standing waves for a NLS on star graphs 
as local minimizers of the constrained energy 

Riccardo Adami* Claudio Cacciapuotij Domenico Fincof Diego Noja 


‘Dipartimento di Scienze Matematiche “G.L. Lagrange”, Politecnico di Torino 
Corso Dnca degli Abruzzi, 24, 10129 Torino, Italy 
e-mail: riccardo.adami@polito.it 

^Dipartimento di Scienza e Alta Tecnologia, Universita dell’Insubria 
Via Valleggio, 11, 22100 Como, Italy 
e-mail: claudio.cacciapuoti@uninsubria.it 

Racolta di Ingegneria, Universita Telematica Internazionale Uninettuno 
Corso Vittorio Emanuele II, 39, 00186 Roma, Italy 
e-mail: d.finco@uninettunouniversity.net 

^Dipartimento di Matematica e Applicazioni, Universita di Milano Bicocca 
Via R. Cozzi, 53, 20125 Milano, Italy 
e-mail: diego.noja@unimib.it 


Abstract 

On a star graph made of A > 3 halflines (edges) we consider a Schrodinger 
equation with a snbcritical power-type nonlinearity and an attractive delta 
interaction located at the vertex. From previous works it is known that 
there exists a family of standing waves, symmetric with respect to the 
exchange of edges, that can be parametrized by the mass (or L^-norm) 
of its elements. Furthermore, if the mass is small enough, then the cor¬ 
responding symmetric standing wave is a ground state and, consequently, 
it is orbitally stable. On the other hand, if the mass is above a threshold 
value, then the system has no ground state. 

Here we prove that orbital stability holds for every value of the mass, 
even if the corresponding symmetric standing wave is not a ground state, 
since it is anyway a local minimizer of the energy among functions with 
the same mass. 

The proof is based on a new technique that allows to restrict the analysis 
to functions made of pieces of soliton, reducing the problem to a finite¬ 
dimensional one. In such a way, we do not need to use direct methods of 
Calculus of Variations, nor linearization procedures. 
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1 Introduction 


The subject of nonlinear dynamics on quantum graphs (or networks, see [HHH 
for an exhaustive introduction) dates back to the seminal papers by Ali Mehmeti 
(see m and references therein), where the focus was on dispersive properties, 
and by von Below who first set variational problems on a network. Since 
then, the interest increased at a growing rate up to the flourishing of results 
of the last decade, motivated by the need for simple models in contexts where 
two features coexist: a basic environment endowed with branches, junctions, 
ramifications, and the presence of non-negligible nonlinear effects. Such models 
range from quantum optics (see, e.g., [19]) to Bose-Einstein condensation (see, 
6-g-) [27] and, for a more comprehensive introduction to physical applications, 
[23] 1. An important part of recent results concerns the study of the Nonlinear 
Schrodinger Equation (NTS), see e.g. [U H] S] 13 [HI HI [1311131 IHl HU 1131 HB] - 
and many of them are concerned with the seek for standing waves, i.e. solutions 
to the NLS that preserve the spatial shape and harmonically oscillate in time 
(namely, solutions of the form $(t) = where is the space profile, 

called stationary or bound state), or even for ground states, i.e., stationary states 
that minimize the NLS energy among all functions with the same L^-norm or 
mass. 

To this regard, in spite of the quasi one-dimensional nature of networks, the 
structure of the family of standing waves, as well as the problem of the existence 
of a ground state, is far richer and more complicated than for the NLS on the 
line (for recent developments in this direction, see [HIH])- 

In the present paper we consider a star graph Q made of N (> 3) halflines 
that meet one another at the unique vertex V (see FiglTj), and study the dynamics 
generated on it by the NLS with an attractive d-interaction of strength —a, 
a > 0, placed at V (for the precise definitions see formulas (I2.1|l . (12.21) 1. It 
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Figure 1: A star graph made of four halfline and a vertex. 

is already known (see mm that for any frequency uj G ,+oo) there 

exists a unique, real stationary state whose associated solution to the NLS 
oscillates at the frequency w and is symmetric under exchange of edges, namely, 
the restriction of to any halfline gives the same function (see Fig|2|). The 
main result of this paper is the orbital stability of and can be expressed as 
follows: 

Theorem 1. On the star graph Q made of N > i halflines intersecting one 
another at the vertex v, consider the Schrodinger Equation (E3D with a focusing 
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Figure 2: A representation of a symmetric stationary state on the star graph 
made of four halflines. 


nonlinearity of power 2/r + 1 and a delta interaction at V with strength —a, with 
0 < p. < 2 and a > 0. 

Then, given ui > a^/N'^, the unique standing wave symmetric under 

exchange of edges, is orbitally stable. 

We recall here that orbital stability is Lyapunov stability for orbits instead 
of states. Indeed, in order to hold for the dynamics generated by the NLS, that 
enjoys phase invariance, Lyapunov stability must be weakened: it cannot hold 
for states, but it may hold for orbits. In other words, Theorem I establishes 
that a solution remains arbitrarily close to the orbit of provided that the 
initial data had been chosen as suitably close to the same orbit. 

With Theorem 1 we complete the analysis carried out in [4], where the 
existence of a ground state at a fixed mass has been established assuming that 
the mass is smaller than a critical value. In that case, orbital stability follows 
from the fact that is a ground state, so the celebrated general result by 
Cazenave and Lions m applies. On the other hand, if its mass is larger than 
the critical threshold, then 4',^ is not a ground state (see point 4. in Section 
ID), therefore one is forced to use a criterion for orbital stability for which it 
is not necessary to assume that is a ground state. In fact. Theorem 3 in 
PP] establishes that a stationary state is orbitally stable if and only if it is a 
local minimum for the energy functional among the functions with the same 
mass. Even though, in general, proving that a function is a local minimum of a 
functional is a difficult task, in our case it is possible to exploit the particular 
structure of the star graph made of N halflines, and Th. 4.1 in [5] which 
states that the minimizer of the NLS energy on the halfline under the mass 
constraint and a nonhomegeneous Dirichlet condition at the origin is given by 
a unique soliton branch (we recall its explicit expression in Sec. P]). This fact 
allows to introduce the so called multi-soliton transformation, that maps almost 
every function on the graph into a function (multi-soliton) made of N pieces 
of solitons, one for each halfline, in such a way that the mass is preserved and 
the energy is lowered. The space of multi-solitons with the same mass, denoted 
by M (see Def. 12.111 . is a finite-dimensional manifold that contains all the 
stationary states. Thus, proving that a stationary state is a local minimum in 
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j\4 requires a finite-dimensional analysis only. Nevertheless, this turns out to 
be not an immediate issue, because the study of the sign of the Hessian of the 
energy requires a certain degree of explicitness. In order to get it, one has to 
further reduce the problem from JV halflines to two, and use the fact that on the 
real line the orbital stability has been already proven by Fukuizumi, Ohta and 
Ozawa in [18) . Finally, owing to the continuity of the multi-soliton transform, 
one shows that the symmetric stationary state is not only a local minimum of 
the energy among the multi-soliton states, but also among all states with the 
same mass. 

The statement of Theorem 1 emphasizes the dynamical content of the result, 
i.e. the orbital stability. It emerges from the preceding discussion that one can 
give a variational version of the same result, that highlights the fact that the 
examined standing wave corresponds in fact to a local minimizer of the energy 
in the appropriate space, which is in general a highly non-trivial goal in the 
Calculus of Variations. For this reason we remark that Theorem 1 can be 
stated also in the following variational version: 

Theorem 1’. (Variational version) On the star graph Q made of N > 3 halflines 
intersecting one another at the vertex v, consider the Schrodinger Eguation (EH) 
with a focusing nonlinearity of power 2p, + 1 and a delta interaction at V with 
strength —a, where 0 < /i < 2 and a > 0. 

Then, given uj > /N"^, the unique positive bound state symmetric under 

exchange of edges, is a strict (up to phase invariance) local minimizer for the 
energy functional associated to the considered evolution equation, constrained to 
the manifold of constant L^-norm. 

In [2 we treated the dual problem, namely the minimization of the action 
functional on the associated Nehari manifold. It was proved that in order to 
have a minimizer, the delta interaction must be strong enough. The link with 
the small mass condition given in [1] can be reconstructed by noting that, given 
a frequency w, the mass of is a monotonically decreasing function of the 
interaction strength a, so that the assumption of large a can be thought of as 
a small mass hypothesis. Through Theorem I one can get rid of every such 
assumptions: 'I'^ is always orbitally stable, irrespective of the mass and of the 
interaction strentgh. Notice that in it was proven that the energy constrained 
at constant mass has also other stationary points than the symmetric states 
above discussed and so the NTS on a star graph admits several families of non 
symmetric standing waves. In [4] it is shown that they are excited states, in the 
sense that their energy is above the energy of the symmetric state. They are 
believed to be saddle points of the constrained energy, which is consistent with 
the fact that the local minimum represented by the symmetric state is above 
the infimum of the energy for large mass. 

Theorem 1 is already known to hold for V = 2 too, i.e. for the NTS with 
an attractive delta interaction on the line, since in that case (which is the 
unique bound state at its frequency and at its mass) is always a ground state 
(see [711I71II1I21]). 

The structure of the paper is the following: in Section [5] we introduce the 
notation, recall the state of the art (in particular we explain the result of Th. 4.1 
in 0 ), give some preliminary results, introduce and describe the multi-soliton 
manifold A4 . In Section [3] we introduce the multi-soliton transformation, that 
allows to reduce the problem to a finite-dimensional setting. In Section 0] we 
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prove that is actually a local minimum for the restriction of the NLS energy 
to j\4. Finally, in Section [5] we extend the minimality property of from the 
multi-soliton manifold to the natural space of the functions with finite energy 
and constant mass, so accomplishing the proof of Theorem 1. 

Acknowledgements. The authors are partially supported by the FIRB 2012 
project “Dispersive dynamics: Fourier Analysis and Variational Methods”, code 
RBFR12MXPO. Furthermore, R.A. is grateful to the PRIN 2012 project “As- 
petti variazionali e perturbativi nei problem! differenziali nonlinear!” and to 
the 2015 GNAMPA project “Propriety spettrali delle equazioni di Schrddinger 
lineari e nonlineari”. C.C. acknowledges the support of the FIR 2013 project 
“Condensed Matter in Mathematical Physics”, code RBFR13WAET. 


2 Setting, notation, previous results 


We follow in notation the previous papers Elia 12 [5]. The basic environment is 
the metric graph Q, defined as the star graph made of N halflines intersecting 
at their origins, where the unique vertex v is located. This convention fixes a 
coordinate system on Q, namely, each edge is identified with the real nonnegative 
halfline [0, +oo) and is endowed with its own coordinate (denoted by xj for the 
j.th edge), while xj = 0 is the coordinate of v, regardless of j. Functions 
dt : Q —> C can be represented as vectors where the component on the j.th 
halfline is a scalar function ipj : [0, +oo) —>■ C. For the sake of clarity, we shall 
often use the column vector representation 




/ -ipi ^ 

i’2 


\ IpN J 


Besides, the component of T on the j.th edge is denoted by ipj or by ('l')j. In 
general, we use uppercase greek letters for functions defined on Q, and lowercase 
greek letters for functions acting on the halfline or on the line. 

We denote by 'k' and T" the vector valued functions with components ijj'j 
and ij)” respectively, the derivatives taken along the coordinate Xj. 

Spaces LP{Q) are naturally defined as the direct sum of N copies of T^(0, +oo), 
while H^{Q) is the direct sum of N copies of H^{0,+oo) with the additional 
condition of continuity at the vertex 


'(/'l(O) = •02(0) = • • • = 0Ar(O) = ^'(v). 


On Q we consider the dynamics generated by the Schrddinger equation with 
power nonlinearity and an attractive point interaction at the vertex, formally 

idt'i'it) = -AT(t)- 

where fi > 0 and the notation |v]/|2/iv[/ is understood in vector representation as 

\ 

/ 




102 P '"02 
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moreover, a > 0 and Sy denotes a delta potential located at the vertex. More 
precisely, 

idt'i’it) = ( 2 . 1 ) 

where H is the linear operator defined as 

N 

D{H) = {d/ e H\g), s.t. V', e i?"(0,+oo), = -avl/(v)} 

1=1 1 1 

(ijvh), = - 

Global well-posedness of (12.11) (for 0 < /r < 2) has been proved first in [I] for 
the cubic case, then in for the case of a more general power, and for the case 
N = 2 (real line) with a general point interaction in [Q. Furthermore, it has 
been also proved that the L^-norm, or mass 

QW = 11^111,2(5) 

and the energy 

m,G) = ^ll^'lli2(5)-^^||iIi||iCi(5)-||di(0)p (2.3) 

are conserved by the flow. 

We shall occasionally make use of the functionals E{-, K+) and E{-, K.), that 
share the formal expression of (ESI) but are evaluated on functions on the halfline 
and on the line, respectively. 

As explained in SecHJ we aim at proving that is a local minimizer for 
E{-,g) with the mass constraint, so, chosen M > 0, our reference space is 

HUG) := H\g)f^{Q{'i>) = M}. 

Let us now recall some preliminary notions, together with some further basic 
definitions and notation. 

1. For every w > 0, the soliton 

4>ui{x) = [{p+ l)a;]^ sech?'(/r-v/wa;) 

is the unique positive square-integrable solution to the stationary NLS equation 

if" (x) + (x) = UJ(p{x). 

As a consequence, the function is a standing wave for the NLS on the 

line with power nonlinearity 2^+1. Besides, for any a; > 0 the mass of the 
soliton is given by 

Mhm = [\l - (2.4) 

M do 

and is a monotonically increasing function of uj. 

2. As proved in [3, for every a > 0, the unique positive stationary solution 
to Eq. (12.11) with frequency uj > a'^ jN"^, symmetric under exchange of edges, is 
given by 

= (jjujiXi + C), i=l,...,N 
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with 


C = —\= arctanh [ —| . (2.5) 

\N^J ^ ' 

Notice that C > 0, so that every halfline hosts a soliton tail, which is a mono- 
tonically decreasing function (see Figl2). 

3. The mass function 


M{uj) := Q(«'^) = f {l-f)^-^dt, (2.6) 

(see formula (5.1) in [1]) is strictly monotonically increasing, and ranges from 0 
(excluded) to +oo as uj goes from o? (excluded) to +oo. Then, in the same 
way as for the soliton on the real line, a; can be interpreted as a relabelling of the 
mass, and for every positive M there exists exactly one symmetric stationary 
state with mass M (see [4]). 

4. As mentioned in SecH] there are some values of M such that the corre¬ 
sponding stationary state is not a ground state. Indeed, in we exhibited a 
sequence G s.t. E{^mQ) —>■ A(())^g,]R) as n goes to infinity (see for¬ 

mula (3.12) in [3), where wr is the unique value of uj such that ||(/)a;|1^2(-R) = M, 
see Eq. (1^ . Roughly speaking, such a sequence is supported on a single edge 
and asymptotically reconstructs a soliton at infinity. Therefore, in order for 
to be a ground state, it must be < E{(j)^,^,R). By explicitly computing 

the involved energies, see also formulas (4.10) and (4.12) in [4], such inequality 
can be rewritten as 

{2-^)wrM < {2 - ^)u;M + afj,{fi +1)'!^ ^uj - . (2.7) 

Furthermore, since ||4',^|||2(g) — M — 11 11 ^2 (R), and from Eqs. (El) and 

(EH), one has the identities 


M 


= ^kLlDl J-i 


(1 - C) 


2 \ 77-1 


dt 


(^ + 1)»‘ i_i 

= N— -2 


(l-t^)M-vt. 




( 2 . 8 ) 


Hence, for M large, the l.h.s. of EH is of order M , as the first term in 

2 

the r.h.s., while the second term in the r.h.s. is of order M , then it can 
be neglected. Expliciting wr and u as functions of M in EH, one has that 
inequality EH amounts to 

{1 — t^)'i^~^dt > — f {1 — f^)'i^~^dt, 

that is violated for M large, as w becomes large too. Then, large mass implies 
that is not a ground state, and, since by Lemma 5.2 in [4], is the 
minimizer of the energy among all stationary states, one concludes that there 
is no ground state if the mass exceeds a critical threshold. 
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The last part of the present section is devoted to the introduction of the 
finite-dimensional manifold to which we shall reduce the problem. 

We preliminarily recall Theorem 4.1 of [5], which establishes that, given 
a,m > 0, there exists a unique couple w > 0, ^ S K such that 

/• + O0 

/ + ^) dx = m, (2.9) 

^0 

Furthermore, the function minimizes the energy 

among the functions in with mass m and whose value at zero equals 

a (Dirichlet constraint), hence, within such class of functions, it is also the 
minimizer of 

By this result one can introduce two functions lu = uj{m,a),^ = ^(rn,a), that 
give the value of the soliton parameters w and ^ such that -|- ^) minimizes 
the functional i?(•,]&+) with both mass and Dirichlet constraints given by Eq. 

(EH). 

As we will show in SecEl the existence of such functions is the starting 
point for the reduction of the problem to a finite-dimensional manifold, whose 
definition is: 


Definition 2.1. Fixed M > 0, we call multi-soliton manifold of mass M, and 
denote by M., the subset of made of functions whose restriction at every 

halfline of Q gives a piece of soliton. 


Remark 2.2. Notice that all functions in A4 are positive. 
Remark 2.3. Every element of A4 has the following form: 

/ </>t.u(Tni.a)(- +C(wi,a)) 

4^Lj{m2,a) (* d” a)) 


<i) 


7n,a — 








where 


( 2 . 10 ) 


• we used the notation m := {mi, m 2 , ■ ■ ■ mN-i), with mi > 0 and < 

M; 

• by the definition of the functions oj and mi is the mass located on the 
bth edge, i.e., 

4>l{m.,a)ix + iimi^a))dx = m^, (2.11) 

and a is the value attained by $^,0 at v, in particular, for all i = 
(1,...,7V-1) 

(^(ar^, a)) = u. (2.12) 

Notice that the global mass constraint M($m,a) = M is guaranteed by the last 
compoment of the vector in the r.h.s. of (12.1011 . 

Since the function depends on N parameters mi ,..., mN-i, a, A4 is a 
iV-dimensional submanifold of H\^{Q). 








Remark 2.4. The real symmetric stationary state 4'^ with (5('kw) = be¬ 
longs to M. Indeed, it is immediately seen that 4'^ = where the vector 

m reads fhi = • • • = ffiN-i = ^ and a is defined as the unique solution to the 
transcendental equation 

~ ^ ^ / Q? \ 

a = K{MlNra){i{^l^^o)) = [{fi + l)uj{M/N, a)] |^1 - J > 

where we observed that ^{M/N,a) = ( and then used (I2.5I1 . 

When dealing with functions in the manifold M, the energy functional be¬ 
comes a real-valued function of N real, positive variables. We emphasize this 
change of point of view by introducing the reduced energy as 

N-l 

Er : {(m, a) G (0,-l-oo)^, s.t. < M} —»■ R 

i=l 

Er{m,a) 

Let us define the function f : [0, M] x R+ —^ R as 
F(m,a) = +^("i,a))||i 2 (R+) 

- +eKa))||^'^t?2(R+) - 

so that the function Ej. decomposes as follows 

N-l N-l 

Er{rn,a) = E{mi,a) + E{M — mj,a). 

i=i j=i 

3 Multi-soliton transformation 

The first step of the proof of Theorem 1 consists in reducing the problem to 
a finite-dimensional manifold, made of functions obtained by gluing together 
pieces of soliton. To this aim, we start by transforming every function in 
that does not vanish at v into a function of AI, in such a way that the value at 
V and the mass at any edge are preserved. 

Owing to Theorem 4.1 in [5] and to the definition of the functions w and 
introduced in SecHJ it is possible to give the following 

Definition 3.1. The soliton transformation Sr] G R^(R+) of a function rj G 
iL^(R+) such that r]{0) ^ 0, is defined as 

Sv ■= +C("l,a)), 

where we denoted m = |? 7 (a;)p dx, a = |r/(0)|. 

By Theorem 4.1 in [8], the function Sg is unique and 


(2.13) 


(2.14) 


(2.15) 


9 






hence, since («S'? 7 )( 0 ) = |? 7 ( 0 )|, Theorem 4.1 in [5] implies 

E{Stj,R+) < E{ri,R+), 

where equality holds if and only if 77 = + ^), for some values oi 9, co, and 

?• 

Furthermore, notice that the soliton transformation acts trivially on pieces 
of soliton, namely + 0 = + 0- 

Finally, the soliton transformation is defined for a > 0 only. However, we 
do not need to cover the case a = 0. The only property we shall use is the 
continuity of S' at the stationary state 4/^. 

Proposition 3.2. The soliton transformation S is continuous from the space 
H\^{Q) to itself. 

Proof. We decompose S in three steps and show continuity at every step. 
Schematically, 

r] (a,m) (w,^) (/)<^(- +^), 

where a = | 7 ?( 0 )|, m = \ri{x)\‘^dx. 

First, Si is continuous because the pointwise value and the L^-norm are 
continuous in i?^(K+). 

As already stated, the fact that the map S 2 is well-defined follows from Theorem 
4.1 in [ 8 ]. The proof that it is also continuous can be made by closely following 
the proof of Theorem 4.1 in [5]. We sketch the procedure since, due to differences 
in the notation, this step could not be straightforward. 

By using the scaling property 



in the identities ( 12.1111 and ( 12 . 121 ) we obtain 

^ _ 2 —pi _ _1^ 

(j)l{x + \/(jj^)dx = moj and = auj . (3-1) 

Putting them together one obtains the identity 

(3.2) 

with 

rOO 

g{z) = {(j)i{z)) ^ / (j)l{x-\-z)dx. 

Jo 

The function g is continuous and strictly monotonically decreasing (for the proof 
of this statement we refer to the proof of Th. 4.1 in [ 8 ] again). We remark that 
the definition of the function (fi in [5] is slightly different from ours, basically 
the two definitions involve different scalings of the hyperbolic secant. This is 
due to the fact that in [5] the authors parametrize solitons through the mass 
instead of the frequency. In spite of that, the argument in Th. 4.1 still applies, 
being based only the asymptotic, monotonicity, and log-concavity properties of 
sech. 

By (13.2p . we infer that the quantity y/iof, is a continuous function of m and 
a. Moreover, for fixed m and a there exists a unique value of y/ujf, such that 
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identity is satisfied. Furthermore, fixed the quantity there exists a 
unique lu such that the second identity in (13.11) is satisfied, and such a w is a 
continuous function of a and m. As a consequence ^ is a continuous function of 
a and m too, and this proves the continuity of 5'2. 

In order to prove the continuity of S 3 , let us fix the couple (w,^) and prove 
that 

lira + C) “ (■ + ^i)llffi(R+) = 0- (3-3) 

Preliminarily, we use the triangular inequality 


0 (• + Ci)IIh1(R+) 

< + ^) — + Cl)llli-i(R+) + + ^ 1 ) — (• + ?l)llfl-i(R+) 

^+00 

< 2||(^a;(* + 0 ^ / 14^10 {X ^ 1 ) — {x + 

JO 

p +00 

+ 2 / \(l}Ux + ii)-(l}L,{^ + ii)\^dx 

JO 

/ + 00 

\4>uj{x) - (l)u^,{x)\'^dx 

-00 


+ 2 



Wuiix) - (t)l,{x)fdx 


(3.4) 


Concerning the first term in the r.h.s., a straightforward computation gives 

0 — (/'wC- + ?l)||i2(R+) < 1C — Cini?^>Llli2(R+)) 

and, repeating the same computation for the derivatives, one gets 

+c) ~ +ci)iifl'i(R+) < 1 C~ ciniC’(.jiiff 2 (R+). (3.5) 

So the first summand in the r.h.s vanishes as Ci approaches C- For the second 
term in the r.h.s. of dH, first observe that, assuming oj /2 < wi < 2 uj, 

- 4>ujiix)\‘^ < 2\4)ujix)\'^ + 2\(j)i^-,^{x)f 

and the last quantity is an integrable function in the variable x. Analogously, 
for the last term in the r.h.s. of inequality (13.41) . one gets 

\JJx)-J^^{xr < 2|</>L(x)p+2|</>U(x)P 


Then, by dominated convergence theorem, one has that the last two terms in 
dH vanish as oji goes to w, thus dISl) is proved and the proof is complete. 

□ 

We can now introduce the multi-soliton trasformation E as the natural gen¬ 
eralization of the soliton transformation S to the star graph Q. 
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Definition 3.3. Given a function $ G H\j{Q) such that <f>(v) ^ 0, the multi- 
soliton transformation of is the function defined on M as 


(!]$), := 

where S is the soliton transformation introduced in Definition AS. 1[ 


Remark 3.4. The multi-soliton transformation S inherits the following prop¬ 
erties from the soliton transformation S: 

1. S is continuous from the space of the functions in H\^(Q) that do not 
vanish at the vertex, to 

2. S preserves the mass distribution on the edges of the star graph, i.e. 
II(^‘^’)jIIl 2(R+) = ll(‘^)illL2(R+), and the absolute value of the function in 
the vertex, namely E$(v) = |$(v)|. 

3. For every $ in the domain of E, 




(3.6) 


where equality holds if and only if $ G Ai, up to a constant phase factor. 

4. The multi-soliton transformation E acts trivially on Al. In particular, 
E4-^ = 4',.. 

4 Local minimality of in M. 

Here we treat the finite-dimensional problem of the local minimality of in 
the manifold A4. In this section we always refer to the reduced energy Ej. (see 
Definition l2.13|) and then, when possible, avoid any reference to functions, that 
can be replaced by points of In particular, according to Remark [521 the 
bound state corresponds to the point P := {m,a) = {M/N,... ,M/N,a). 
However, in order to prove even this finite-dimensional minimality, along the 
proof we must come back to function representation and make use of a well- 
known result of the Grillakis-Shatah-Strauss theory on stability of standing 
waves (see Theorem 3.4 in [5D]), that we straightforwardly apply in order to get 
inequality 

Proposition 4.1. Fixed M > 0, the point P is a strict local minimum for the 
function Er{m,a) defined in (12.131) . 

Proof. First, notice that P is an internal point for the domain of Er, therefore it 
suffices to show that P is a stationary point for E^ and that the Hessian matrix 
of Er evaluated at P is positive definite. ^ 

The fact that 4/^ is a stationary point immediately gives that P is a sta¬ 
tionary point for Er, so we turn to the study of the sign of the Hessian matrix. 
By straightforward calculations. 



(4.1) 
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with 1 < i, j < N — 1. ^ 

Consequently, the Hessian matrix of Er computed at P is a N x N block matrix. 
The high {N — 1) x (TV — 1) left block, that we call Hi, reads 


d^F fM 


111 := 


dm? \ N ’ 


(J + I) 


where J is the matrix with all elements equal to one. By elementary linear 
algebra, one immediatly finds that Hi has two eigenvalues: 
multiplicity 1, and multiplicity N — 2. Therefore, Hi is positive 

definite if and only if 


d'^F 

dm^ \ TV ’ ) 


> 0 . 


(4.2) 


The second diagonal block, denoted by H 2 , is 1 x 1 and reduces to (^, a), 

so it is positive defined if and only if 


d^F (M \ 
da? \ TV ’ y 


> 0 . 


(4.3) 


Thus, in order to prove that P is actually a local minimum of the function Er, 
one has to prove 621) and 6S1)- 

To prove (14.21) . consider the curve in the manifold M given by ^rn{t),aj with 

, , M , , , . M , , 

mi(t) = — + t, m2[t) = ■■■ = mN-i{t) = ^ ^ (“£:£) 

and define the function 

f{t) := Erim{t),a). 

By definition of partial derivative and using iU, one has 


no) = ?FnP) = 2^ (“ si 

' dm? ' )' 

So condition (14.21) reduces to /"(O) > 0. Now, let us write the curve in terms 
of elements of Af: 




^ 4>Lj{M/N+t,a){.' P + \ 

4>uj(M/N,a) {' + C(-^/-^i«)) 

\ (l>ui(M/N-t,a){- + HPI/N — t,a)) J 


It transpires that ^m{t),a varies with t in the first and the N.th components 
only. Since the total mass is conserved, as t changes there is a transfer of mass 
from one edge to the other without changing the value at the vertex. We define 
the operator 

r : L^{g) L‘^{R) 

acting as rS = rj, with 
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where, as usual, we denoted by (E!)j the component of the wave function S on 
the j.th edge. In other words, 77 is the function on the line obtained from S by 
matching together edges 1 and N and neglecting all the others. Therefore, 


f{t)-EriP) = Erim{t),a) - EriP) 

F{m,{t),a)+F -iVF ("^,0 


= E2a/NiT^rn{t},ay^) “ ^2ct/A7 ; I^) ; 


(4.4) 


where we exploited the definition (I2.14|) of the function F and introduced the 
functional 

E2 o./n{u,R) := Iwu'Whm - - ;^l«(0)l' 

acting on Notice that £' 2 a/Ar(',R) is the functional representing the 

energy associated to the focusing Schrodinger equation with power nonlinearity 
2/r+ 1 and a delta interaction placed at the origin with strength 2a/N (see [TH] 
and [6]). 

Furthermore, 

(T4'^)(a:) = [(7r + l)a;]37rcosh"'^ (7X\/w(|a;| +C)) , 

with C given by Eq. (E^, which is the ground state of the NTS on the line 
with a delta interaction located at the origin, of strength a' = ^. In that case 
(see [H] and m): such a state is known to be a stable global minimizer of the 
constrained problem, so it falls into the scope of Theorem 3.4 of m- In the 
notation of that theorem, T(s{u)) = 1 since we are dealing with real functions. 
Therefore, for e sufficiently small, the theorem yields 

f(t)-Er(P) = E2a/w(T$™(t),S,K) - E2a/Jv(T^'a;,R) 

> c||T$„(t),s-r4'^||i2(R), (4.5) 

for any t G (—e,e). We claim that there exists a positive constant c such that, 
for e small enough, 

- T4'a;||i2(R) > C (4.6) 

This gives the bound f{t) — Er{P) > ct^, which, together with f{0) — Er{P) = 0 
and /'(O) = 0, implies /"(O) > c > 0, and concludes the proof of (14.21) . 

It remains to prove the claim dMl)- To this aim, we write 

\\T^m(t),a — T^cj|li2(R) 

{(l)uj{M/N+t,a){x + ^{M/N + t, a)) - (j)uj{M/N,a){x + ^{M / N ,a))Ydx 

.'OO 

+ / {(l)u:{M/N-t,a)ix + ^{M/N - t,d)) - (l)^(^M/N,a)ix + ^iM/N,d))ydx 

Jo 
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and note that 


||t$. 


m(t),o 


— 


“llz,= (R)lt =0 


= 0 and 

at 


m(t),a 


- 


Il2(R) 


= 0 . 


t=o 


To conclude the proof of (ITO it is enough to show that 




> c. 


i =0 


We compute the second derivative at t = 0 and obtain (we omit the dependence 
of uj and ^ on m and a) 




— T'h, 




t =0 


= 4 


dui dm dx dm ' 


(4.7) 


(m,a) = (M/N,a) 


Next we prove that the integrand is not identically equal to zero, which in turn 
implies that the second derivative in t = 0 is strictly positive. 

In Eq. (13.21) we set z = ^/cJ^ and take the derivative with respect to m, this 
gives the identity 

dff dz 1 
dz dm ’ 

which tells us that ^ 0. Since 

dz ^ duj ^ d^ 

dm 2y/uj dm ^ dm ’ 

we conclude that ^ and cannot be both equal to zero (recall that ^ 7 ^ 0 
whenever a 7 ^ 0). Since the functions ^^(x) and ^q^(x) are linearly indepen¬ 
dent, the integrand in Eq. (03) does not vanish identically. 

In order to prove 631) one proceeds analogously. First define a curve ^fh,a{t) 
in the manifold AI, with 


a(t)=a-\-t, t£{—e,e). 

Then, introduced the function h(t) := Er(m,a{t)), one immediately has 


h"{0) 


d'^Er 

da? 


(P) 


N 


d‘^F 

da^ 



Hence, to prove (14.31) it is enough to show that h"(Q) > 0. Notice that, for any 
t, the function is changing symmetrically on any edge, so that, using 

Theorem 3.4 of m again, 


h(t)-Er(P) = ArF(^^,a(t)^ 

= -^(-E'2a/Ar(T‘f’m,a(t),R) ” ^'2a/iV (t'^w , R), ) 
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To prove the inequality h''{0) > 0 it is enough to show that ||T$^ ^(t) ~ 
> ct^ for t small enough. Arguing as in the proof of (14.61) . we 

conclude that this is certainly true if ^ and ^ are not both equal to zero. 
To see that this is actually the case, we take the derivative of Eq. dm with 
respect to a and obtain an identity which is not compatible with = 0. This 
in turns implies that or 0 or || 0. 

□ 

Remark 4.2. As an alternative formulation of Proposition 14.11 we have that 
'I'tj is a strict local minimizer for the restriction of E{-,Q) to the manifold Ai. 


5 Orbital stability 

The next step consists in passing from the local minimality of in Ai to the 
local minimality of 4'^^ on This step is immediate once one considers 

that: first, the reduction from H]^{Q) to Ai through the transformation S 
lowers the energy level; second, 4',^ is invariant under E; third, E is continuous 
at 4^^. Some care must be dedicated to the fact that multiplying by a constant 
phase factor does not lower the energy. Anyway, one has 

Proposition 5.1. 4t,^ is a strict (up to multiplication by phase) local minimizer 

ofE{;g) inHl,[g). 

Proof. First we prove that 4t,^ is a strict local minimizer among real functions 
in H\j{g). According to Remark 14.21 4/,^ is a strict local minimum of E{-,g) 
restricted to Ai. This means that there exists e > 0 such that, if 4> y^ 4>,^ is a 
real element of H\j{g) with 4>(v) y^ 0, and ||E$ — ^t.j\\m{g) < £, then 

Ei-^^,g)<E(E<i>,g), ( 5 . 1 ) 

where equality holds if and only if E4> = 4',^. 

Moreover, by Rermark 13.41 we know that E is continuous at 4>(^, then there 
exists 5 > 0 such that if 0 < ||4> — 4'^||//i(0) < 6 , then 

e > ||E4) — E4'(^||^i(g) = ||E4) — 

where we used the invariance of 4t,^ under the action of E. But then, using 
inequalities (EU and dSH), 

A(4'^,a) < A(E$,a) < A($,a), 

where the first inequality becomes an equality if and only if E$ = 4'(^, while 
the second one becomes an equality if and only if $ = E4>. Then, we proved 
that there exists a (5 > 0 such that if 0 < ||4) — 4',^||/yi(0) < S, then A(4>, g) < 
E{'i’,g), so we have that 4>,^ strictly minimizes E{-,g) locally among the real 
functions in E[\^{g). Of course, extending the analysis to non-real functions, 
cannot be a strict minimizer due to phase invariance of the energy functional: 
E{e^^'^^,g) = E{'i!i^,g). However, for any 4> outside the phase orbit of 
whose distance from the orbit is less than 5, one has E(^,g) > A(4',^,C/). 
Indeed, there exists 9 G [0,27r) s.t. 

<5 > ||4> - = ||e-*«4> - 4-^||yyi^(e) > |||e-*^4>| - 


16 


so that, since |e is real and different from 


< Ei\e-^^<i>ig). (5.2) 

On the other hand, a straigthforward computation gives 
E{\e-^^^\,g) < £;(e-*®$,e), 
that, together with (ESI), concludes the proof. 

□ 


So we can conclude by invoking Grillakis-Shatah-Strauss theory. 

Proof of Theorem 1. Owing to Theorem 3 in |20) . the local minimality property 
in E[\j{g) is equivalent to orbital stability. The proof is complete. □ 
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